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ABSTRACT 



Context. Ultracompact (UC) H II regions constitute one of the earliest phases in the formation of a massive star and are characterized 
by extreme physical conditions (Go>10'' Habing field and n>10* cm"^). The UC HII Mon R2 is the closest one and therefore an 
«"^ ' excellent target to study the chemistry in these complex regions. 

'"^ Aims. Our goal is to investigate the chemistry of the molecular gas around the UC H II and the possible variations due to the different 

local physical conditions. 

(~* • Mefhods. We carried out a 3mm and 1mm spectral survey using the IRAM 30-m telescope towards three positions that represent 

O ' different physical environments in Mon R2: (i) the ionization front (IF) at (0",0"); two peaks in the molecular cloud (ii) MPl at the 

I , offset (+15",- 15") and (iii) MP2 at the farther offset (0",40"). In addition, we carried out extensive modeling to explain the chemical 

O ■ differences between the three observed regions. 

Results. We detected more than thirty different species (including isotopologues and deuterated compounds). In particular, we de- 
C/3 ' tected SO^ and C4H confirming that UV radiation plays an important role in the molecular chemistry of this region. In agreement 

^ I with this interpretation we detected the typical PDR molecules CN, HCN, HCO, C2H, and C-C3H2. There are chemical differences 

between the observed positions. While the IF and the MPl have a chemistry similar to that found in high UV field and dense PDRs 
like the Orion Bar, the MP2 is more similar to lower UV/density PDRs like the Horsehead nebula. Our chemical modeling supports 
this interpretation. 

In addition to the PDR-like species, we also detected complex molecules such as CH3CN, H2CO, HC3N, CH3OH or CH3C2H that 
are not usually found in PDRs. A wealth of sulfur compounds, CS, HCS^, C2S, HtCS, SO and SO2 and the deuterated species DCN 
and C2D were also identified. The origin of the complex species requires further study. The observed deuteration fractionations, 
[DCN]/[HCN]~0.03 and [C2D]/[C2H]~0.05, are among the highest in wami regions. 

Conclusions. Our results show that the high UV/dense PDRs present a different chemistry from that of the low UV case. Some 
CO ■ abundance ratios like [CO^]/[HCO+] or [HCO]/[HCO+] are good diagnostics to differentiate between them. In Mon R2 we have the 

C,^ ■ two classes of PDRs, a high UV PDR towards the IF and the adjacent molecular bar and a low-UV PDR which extends towards the 

OQ ' north-west following the border of the cloud. 
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;_! 1. Introduction exist and that lifetimes should be one to two orders of magni- 
Cd tude larger Several models have been proposed to explain this 
During the earhest stages of star formation, the UV radiation "lifetime paradox" but all have shortcomings. The paradox could 
from the newly -born star ionizes the most exposed layers of the be resolved if the molecular gas in which an O star forms is 
host molecular cloud, creating a layer of ionized gas (constituted denser and warmer than previously beUeved, resulting in an ini- 
mainlybyHii) and originating a so-called photo-dissociation re- tial Stromgren sphere much smaller than originally estimated, 
gion (PDR). This suggestion is supported by observations which show that in 
Ultracompact HII regions (UC Hit) are defined as regions of dense molecular cloud cores densities of ~ 1 0^ cm'^ and temper- 
ionized gas with diameters smaller than ~0. Ipc (see Churchwell ^^ures of ~ 100 K are not atypical (Rizzo et al. 2005). The study 
2002 for a review). UC Hii regions are expected to expand at of the physics and chemistry of UC Hii and of the surrounding 
velocities of the order of the sound speed (10 km s') until PDRs is the key to understand this evolutionary stage of a mas- 
reaching equilibrium at dimensions of a few pc. In regions with ^^^^ ^'^^■ 

density of n ~10^ cm"^, Hii regions should remain ultracom- UC Hii are characterized by extreme UV radiation field 

pact for ~3000 yr and only a few dozen should exist in the (G()> 10^ in units of the Habing field) and gas densities often 

Galaxy. However, observations suggest that many more UC Hiis higher than > 10* cm"^. UC Hii may represent the best exam- 
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pie of highly UV-irradiated PDRs and can be used as a template 
for other more complex systems. Highly UV-irradiated PDRs are 
found in astrophysical environments of high interest such as the 
surface of circumstellar disks and in the nuclei of starburst galax- 
ies. Unfortunately, our knowledge of the chemistry and physics 
of UC Hii is still far from complete. 

Monoceros R2 (Mon R2) hosts the closest (d=830 pc; Herbst 
& Racine 1976) and brightest galactic UC Hii, and the only one 
that can be spatially resolved in the mm domain with single- 
dish telescopes. It is therefore the best case-study to determine 
the physical and chemical evolution of massive star forming re- 
gions. The central, spherical UC Hii was created by the inter- 
action of the central infrared source IRS 1 (Wood & Churchwell 
1989) with its host molecular cloud. The formation of the BOV 
star associated to IRSl created a huge bipolar outflow (~15' - 

3.6 pc long, Massi, Felli & Simon 1985; Henning, Chini & Pfau 
1992, Tafalla et al. 1994). Later CO 3^2 mapping of the re- 
gion showed that this very extended outflow is now inactive. 
The highest velocity gas is detected towards IRS 3 suggesting 
that this young star is associated with a compact (<14") bipolar 
outflow (Giannakopoulou et al. 1997). 

Further studies at mm-wavelengths (Giannakopoulou et al. 
1997; Tafalla et al. 1997; Choi et al. 2000; Rizzo et al. 2003, 
2005) showed that the UC Hii region is located inside a cavity 
bound by a dense molecular ridge at the south-east. The peak of 
this molecular ridge is located at the off'set (-1-10", -10") relative 
to the peak of the ionized gas, and presents a molecular hydro- 
gen column density of 2-6x10^^ cm"^. The detection and subse- 
quent analysis of the millimeter lines of the cyclic hydrocarbon 
species C-C3H2 with the IRAM 30m telescope revealed gas den- 
sities of few 10*" cm"^ (Rizzo et al. 2005) in the PDR surround- 
ing the Hii region. Recently, Spitzer observations of the H2 ro- 
tational lines and the polycyclic aromatic hydrocarbons (PAHs) 
mid-infrared bands have revealed a thin layer («=4xl0^ cm"-', 
A^(//2)=lxl02' cm-2) of warm gas (rk=574(+20) K) between 
the ionized gas and the dense molecular gas traced by previous 
millimeter observations (Berne et al. 2009). They interpreted the 
differences between the spatial distribution of the emission of 
the H2 rotational lines and the emission from PAHs around the 
Hi: region as the consequence of variations of the local UV field 
and density in this molecular gas layer. 

In this paper, we present a mm molecular survey towards the 
UC Hii region Mon R2. Our goal is to investigate the molecular 
chemistry in the PDR that surrounds the UC Hii region and try 
to construct a pattern to interpret future observations. 

2. Observations 

Figure 1 displays a summary of the morphology of Mon R2 at 
both mm ('^CO 2— >7) and mid-IR wavelengths. As expected, 
the maximum of the emission at 8 //m is found towards the Hii 
region and the infrared star IRS 3. Extended emission is detected 
around the Hii region that is associated with the PAH band at 

7.7 jum. This emission does not distribute uniformly around the 
Hi! region. It extends more than one arcminute away from IRS 1 
towards the north-west whereas it declines very rapidly towards 
the south-east. 

We performed a 3mm and 1mm spectral survey towards 
three selected positions that are expected to represent dif- 
ferent physical and chemical conditions. Our center position 
(RA=06h07m46.2s, DEC=-06°23'08.3" (J2000)) corresponds 
to the position of the ionization front (IF) and that of the ioniz- 
ing star IRS 1. This off'set together with the (-h15", -15") form 
a strip across the molecular bar with a sampling of «; 1 beam (see 



Table 1. Summary of 30m observations. 



Offset 


Freq 


HPBW 


rjMB 


rms' 


t 




(GHz) 


(") 




(mK) 


(min) 


IF 


[83.746, 87.461] 


29 


0.81 


6 


55 


(0",0") 


[87.517, 91.232] 


29 


0.80 


4 


46 




[91.227,94.942] 


29 


0.80 


4 


46 




[103.197, 106.913] 


24 


0.80 


3 


46 




[106.907, 110.622] 


24 


0.80 


4 


46 




[204.583, 208.294] 


12 


0.63 


20 


41 




[207.184,210.900] 


12 


0.63 


13 


74 




[216.344.220.059] 


11 


0.63 


15 


28 


MPl 


[83.746, 87.461] 


29 


0.81 


3 


55 


(15",-15") 


[87.517, 91.232] 


29 


0.80 


3 


55 




[91.227,94.942] 


29 


0.80 


4 


55 




[103.197, 106.913] 


24 


0.80 


3 


55 




[106.907, 110.622] 


24 


0.80 


3 


55 




[204.583, 208.294] 


12 


0.63 


20 


41 




[207.184,210.900] 


12 


0.63 


13 


64 




[216.344, 220.059] 


11 


0.63 


10 


28 


MP2 


[83.746, 87.461] 


29 


0.81 


4 


111 


(0",+40") 


[87.517,91.232] 


29 


0.80 


3 


46 




[91.227,94.942] 


29 


0.80 


3 


46 




[103.197, 106.913] 


24 


0.80 


3 


46 




[106.907, 110.622] 


24 


0.80 


4 


46 




[204.583, 208.294] 


12 


0.63 


20 


27 




[207.184,210.900] 


12 


0.63 


20 


24 




[216.344.220.059] 


11 


0.63 


10 


55 



rms in units of T^s with WILMA (spectral resolution of 2 MHz). 



Fig. 1). For simplicity, we will refer to this offset as MPl here- 
after. The farther offset (0",+40") (hereafter MP2) is within the 
extended PAH emission towards the north-west. 

The observations were carried out with the IRAM 30m tele- 
scope at Pico Veleta (Spain) during July 2009. To optimize the 
observing time, we used the dual sideband offered by EMIR. 
This way we can simultaneously observe 8 GHz in the lower 
sideband and 8 GHz in the upper sideband, separated by a gap 
of 8 GHz. The rejection from one band to the other one has 
been measured by the IRAM staff to be *13 dB. This observing 
mode allows to cover a band of 16 GHz at 3mm with only one 
tuning. Initially we observed 16 GHz single polarization (from 
87.517 to 94.942 GHz and from 103.197 to 110.622 GHz) us- 
ing the dual sideband of the 3mm receiver. Afterwards we com- 
pleted our survey by observing simultaneously the 3mm and 
1mm bands (4 GHz at 3mm and 4 GHz at 1mm) and integrat- 
ing the two polarizations. The 3mm and 1mm EMIR receivers 
were centered at 85.339 GHz and 217.937 GHz respectively. In 
the 1 mm receiver, we need to drop one of the polarizations due to 
the poorer quality of the baselines. Two additional settings were 
used to observe the SO+ 9/2^7/2 (207.800 GHz) and SO 54^52 
(206.176 GHz) lines. Simultaneously with the observations of 
these lines at 1mm, we used the 3mm receiver to obtain the high 
spectral resolution spectra of the HCO+ 7->0, HCN 1^0, C2H 
l—>0 and C4H 11—>10 lines shown in Fig. 2. 

The observed spectral ranges, telescope beam efficiencies, 
half power beam width (HPBW) and rms of our observations are 
listed in Table 1 . The observing procedure was position switch- 
ing with the reference located at an offset (4-400", -400"). To 
cover these spectral ranges we used the WILMA autocorrelator 
which provides a spectral resolution of 2 MHz (5.4-6.8 km s"' 
in the 3mm band). This spectral resolution does not allow to re- 
solve the lines at 3mm and therefore it gives information only on 
the integrated line intensity. The WILMA spectra and the line 
integrated intensities are shown in Appendix A. Selected lines 
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Fig. 1. Left: Map of the integrated intensity between -5 and 15 km s ' of the '^CO 2— >7 line observed at the IRAM 30m telescope 
(lower contour at 90 K km s"' and steps of 10 K km s"' , from Pilleri et al. 2012a). The three obsei-ved positions, IF, MPl and MP2, 
are indicated. The beam of the 30m at 3mm towards the 3 observed positions is also drawn. Right: IRAC 8 yum map towards Mon2. 
The extended emission is produced by the PAH bands at 7.7 //m. This figure shows the existence of an extended PDR in the cloud. 



were observed simultaneously with the VESPA autocorrelator, 
providing a spectral resolution of around 40 kHz. These high 
spectral resolution spectra are shown in Fig. 2 and Gaussian fit 
parameters to SO, SO2, SO^ and C4H lines are given in Table 3. 
All the intensities are in unit of main beam brightness tempera- 
ture (Tmb). 



3. Results 

3.1. Line identification 

To identify the lines, we have used three different catalogues, the 
JPL (http://spec.jpl.nasa.gov), the CDMS (http://www.astro.uni- 
koeln.de/cdms/catalog) and that from J. Cernicharo (private 
communication) The observations were performed with 2SB re- 
ceivers with an image band rejection of ~13dB. Therefore, we 
have some bright lines from the image band in our spectra (see 
Fig. Al, A2, A3). In Appendix A, we show the spectra of our 
complete mm survey, and give a list of the detected lines (>3-cr), 
their identification and line integrated intensities (Table Al). We 
also list the >4-cr unidentifed lines and the corresponding fre- 
quencies in the signal and image bands. For clarity, upper lim- 
its are not included in Table Al. Since the molecular lines are 
unresolved, the upper limit to the integrated intensity is simply 
3 X nns x channel velocity width. The Gaussian fit parameters to 
the recombination lines are given in Table A2. 

We have identified 87 lines, out of a total of 105 lines, to- 
wards the (0",0") position. The large number of recombination 
lines shows the presence of ionized gas. In addition, we have 
also a wealth of lines of complex molecules, more typical of 
warm and dense molecular clouds. Towards the MPl, we have 
detected 101 lines, and identified all but 5 lines. At this posi- 
tion the recombination lines are very weak. Towards the MP2 
we have detected a total of 99 lines with only 8 unidentified 
lines. We have tentatively detected the H2CN 5o,2^2o,2 line to- 
wards the MPl and MP2. H2CN was firstly detected in TMCl 
by Ohishi et al. (1994). If confirmed, this region would be the 
second one in which this radical has been detected and would 



corroborate that Mon R2 is a source with an extremely rich car- 
bon chemistry. 

In Table 2, we show a summary of the molecular species 
detected in our survey towards the three positions. We distin- 
guish between certain and tentative detections. Tentative detec- 
tions correspond to species with only one line detected and at a 
level of 3-4-cr. In these cases, we checked that there were not 
any other intense (and undetected) line of the same species lying 
in the spectral ranges covered by our survey. In Table 2 we added 
the reactive ions CO^ and HOC^ previously detected by Rizzo 
et al. (2003) because of their importance to interpret the PDR 
chemistry. Rizzo et al. (2003) did not detect these ions towards 
the MPl and they did not search for them in the MP2. Like CO+ 
and HOC^, the reactive ion SO^ is expected to present very low 
fractional abundance in the shielded part of the cloud. In this 
work, SO^ has been searched in the three positions but only de- 
tected towards the MP. Many molecules have been detected in 
the three positions. Some are well-known tracers of PDRs like 
HCN, HCO, C2H, C-C3H2 and C4H. Others, like CH3CN and 
CH3OH, are complex molecules usually not found in PDRs. At 
the end of Table 2, we show those species with certain detections 
only towai-ds the MPl and the MP2 peaks, such as C4H, HCO, 
SiO, C2D and SO2. The sulphuratted carbon-chain C2S has been 
detected in the MP2 and tentatively in the IF. HNCO has only 
been detected in the MP. 

The detection of recombination lines, PDR tracers and com- 
plex molecules is consistent with the interpretation of the PDR 
as an expanding envelope around the UC Hii region firstly pro- 
posed by Rizzo et al. (2005) and recently confirmed by Fuente 
et al. (2010) and Pilleri et al. (2012a). Even in this situation, we 
expect to see some chemical differences because of the different 
incident UV fields and mixing ratios of the different PDR layers. 



3.2. Higin resolution observations: Detection of SO'' 
CaH 



and 



In Fig. 2 we show the spectra of the HCO+ 7^0, HCN 7 -»0, SO 
54^i2, SO+ 9/2->7f2, C2H 7^0, and C4H 77^70 lines. This is 
the first detection of SO^ and C4H towards Mon R2. These lines 



Tmb(K) 
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Fig. 2. High resolution spectra observed with the 30m telescope towards the positions (0",0"), (+15", -15") and (0",+40"). All 
the lines were observed with a spectral resolution of ~40 kHz. The spectra of the SO^ 9/2^>7/2 and C4H 11—>10 lines have been 
rebinned to obtain a higher signal to noise ratio. 



were observed with the ~40 kHz spectral resolution provided by 
the VESPA autorrelator Gaussian fits to the SO, S0+ and C4H 
lines are shown in Table 3. 

There are significant differences among the profiles of the 
lines shown in Fig. 2 which testify the existence of several gas 
components. A clear absorption is seen in the HCO^ 7— >0, HCN 
7— >0 and C2H l—>0 lines towards the IF position. This self- 
absorbed feature lies at a velocity of ~10.8 km s ' which is the 
central velocity of the SO 54— >43 line emission. This suggests 
that the SO 54^43 emission is coming from the cloud around 
the Hii region. Towards the MPl the profiles of the SO 54—^43 
and SO^ 9/2— > 7/2 lines are wider and centered at ~10.0 km s"^ 
The line profiles towards the MP2 are even more complicated 
(see Fig. 2). Similarly to what happens in the other positions, 
there is a self-absorbed feature at ~10.8 km s"' which is the ve- 
locity of the external cloud, but the HCN 1 — >0 line also presents 
self-absorption at ~9 km s '. The C2H 1-^0 and C4H 11^10 
lines show wide profiles, with a red-shifted wing up to velocities 
of ~18 km s '. The similarity between the profiles of the C2H 
l—>0 and C4H 11—>10 lines supports the interpretation of these 
two species share a similar chemistry and come from the same 
region. Since C4H is found to be abundant in low UV PDRs 



and has not been detected in the high velocity wings of bipolar 
outflows (see e.g. Bachiller & Perez-Gutierrez 1997), these high 
velocity wings are very likely tracing the PDR around the Hii 
region (Fuente et al. 2010, Pilleri et al. 2012a). The existence of 
several bipolar outflows in this regions makes difficult, however, 
to disentangle the origin of this high velocity gas. The profile 
of the SO 54— >43 line is, however, Gaussian and centered at the 
systemic velocity suggesting that its emission originates in the 
dense cloud. 

4. Gas Physical Conditions and column densities 



To derive accurate column densities it is important to have a 
good estimate of the gas physical conditions. For the species 
Ci**0, SO, i^CN, SiO, C2S, SO2, H2CO, H2CS, C-C3H2, HC3N, 
C4H, CH3OH, CH3CN and CH3C2H we detected several tran- 
sitions and a multitransitional study was possible. We used the 
rotational diagram technique to derive beam averaged rotation 
temperatures and column densities. This technique is based on 
the assumption of optically thin emission for all the lines of a 
given species. For optically thin emission, the integrated inten- 
sity of each line is proportional to the population of the upper 
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Table 3. SO, SO2, SO+,C4H observations. 



Line 




HPBW(") 


Freq(GHz) 


I(Kkms-') 


v/„(kms ') 


Av(km s ' ) 


Tms(K) 


Offset (0",0") 
















SO 


54^43 


12 


206.176 


0.70(0.05) 


10.56(0.08) 


2.55(0.27) 


0.26 


SO+ 


912-^712 (f) 


12 


208.965 


0.11(0.05)' 


10.09(1.00) 


3.87(1.72) 


0.026 


SO+ 


912^712 (e) 


12 


208.590 


<0.16 








S02 


52,2^2u 


12 


208.700 


<0.16 








C4H 


11-^10 


24 


104.666 


0.09(0.02) 


9.77(0.46) 


3.56(1.11) 


0.025 


C4H 


11^10 


24 


104.705 


0.11(0.02) 


10.17(0.39) 


3.77(0.80) 


0.028 


Offset (+15" -15' 


') 














SO 


54^43 


12 


206.176 


1.16(0.04) 


10.13(0.04) 


2.59(0.09) 


0.42 


SO+ 


9/2-^7/2 (f) 


12 


208.965 


0.26(0.03) 


10.17(0.12) 


2.06(0.21) 


0.12 


SO+ 


9/2^7/2 (e) 


12 


208.590 


0.13(0.04) 


10.20(0.21) 


1.47(0.53) 


0.08 


S02 


^2,2— >2l,l 


12 


208.700 


<0.07 








C4H 


11^10 


24 


104.666 


0.09(0.03) 


8.70(0.81) 


5.14(1.64) 


0.017 


C4H 


11^10 


24 


104.705 


0.14(0.03) 


9.78(0.63) 


5.90(1.17) 


0.022 


Offset (0", +40") 
















SO 


54^4, 


12 


206.176 


1.19(0.03) 


9.65(0.02) 


1.70(0.06) 


0.66 


S0+ 


9/2^7/2 (f) 


12 


208.965 


<0.17 








S0+ 


9/2^7/2 (e) 


12 


208.590 


<0.17 








S02 


^2,2 ^2 1,1 


12 


208.700 


<0.22 








C4H 


11^10 


24 


104.666 


0.14(0.02) 


6.53(0.07) 


1.11(0.16) 


0.12 


C4H 


11^10 


24 


104.705 


0.09(0.02) 


7.99(0.09) 


0.90(0.18) 


0.10 



' Tentative, <3-cr 



Table 2. Detected molecules. 



Table 4. CoUisional coefficients. 



IF(0' 


',0") 


MPl (+15",-15") 


MP2 (0",+40") 


Detected 


Tentative 


Det. 


Tent. 


Det. Tent. 


CO+ 










HOC+ 




SO+ 






"CO 




"CO 




"CO 


C"*0 




c"o 




c"*o 


SO 




SO 

"CN 
C'^N 




SO 

"CN 
C"N 


"CS 




"cs 




"cs 


N2H+ 




N2H+ 




N2H+ 


HCS+ 




HCS+ 




HCS+ 


HNC 




HNC 




HNC 


HCN 




HCN 




HCN 


H"CN 




H"CN 




H"CN 


HN"C 




HN'^C 




HN"C 


HC'N 




HC'^N 




HC'^N 


DCN 




DCN 




DCN 


HCO+ 




HCO+ 




HCO+ 


H"CO+ 




H"CO+ 




H"CO+ 


C2H 




C2H 




C2H 


C-C3H2 




C-C3H2 




C-C3H2 


H2CO 




H2CO 




H2CO 


HC3N 




HC3N 




HC3N 


H2CS 




H2CS 




H2CS 


CH3CN 




CH3CN 




CH3CN 


CH3OH 




CH3OH 




CH3OH 


CH3C2H 




CH3C2H 




CH3C2H 




C4H 


C4H 




C4H 




C2S 


HCO 

SiO 

C2D 

SO2 

HNCO 


H"NC 
H2CN 
NH2D 


HCO 

SiO 

C2D 

SO2 

C2S 

H2CN 



c"*o 


Yang et al. 2010 


SO 


Lique et al. 2006 


cs 


Lique et al. 2007 


SiO 


Dayou & Balanca, 2006 


HCN, HNC 


Sarrasin et al. 2010 


HCO+ 


Flower, D.R. 1999 


HCS+ 


Monteiro 1984 


N2H+ 
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level of the corresponding transition. In this case, a rotational di- 
agram provides good estimates of the rotation temperature and 
total column density (see e.g. Schloerb et al. 1983). In the case 
of optically thick lines, the integrated intensity of the line is not 
proportional to the population of the upper level anymore and 
the rotation temperature cannot be derived with this technique. 

The lines of the most abundant isotopologues are optically 
thick and in many cases the lines are self-absorbed (see Fig. 2). 
For this reason, we only made rotational diagrams for the rarer 
isotopologues C^^^O, H'^CO+ and H'^CN. The column densi- 
ties of CO, HCO^ and HCN were derived assuming the isotopic 
ratios '^C/'^C=50 (Savage et al. 2002 and references therein) 
and '^O/'^O=500. The information of our survey was completed 
with the H'3CO+ 5-^2, H^^CN 3^2 and C2H 3^2 lines ob- 
tained from our unpublished HERA maps (Pilleri et al. 2012b, 
in preparation) after degrading their angular resolution to that of 
the observed 3mm transitions. The resulting integrated line in- 
tensities are shown at the end of Tables Al. In these cases, the 
rotational temperatures and column densites derived from the 
rotational diagrams are not affected by the unknown beam fill- 
ing factor. For the other cases, we assumed a beam filling factor 
of unity for the 3mm and 1mm lines. This is reasonable taking 
into account that the molecular emission extends over a region 
of about 2'x2' (see Fig. 1), which is larger than the beam of the 
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30m telescope at both 3mm and 1mm wavelengths but it would 
overestimate the rotational temperature if the size of the emis- 
sion is smaller than 29". Comparing the rotational temperatures 
thus derived with those from H'^CO+, H'^CN and C2H lines we 
estimate that the uncertainty in the column density estimates be- 
cause of the unknown filling factor is of a factor of 3. 

Where possible we estimated the molecular column densi- 
tites and the molecular hydrogen density using the large velocity 
gradient radiative transfer code by Cernicharo et al. (2006). The 
collisional coefficients adopted in our calculations are shown in 
Table 4. For '^CS, Hi3CO+, H^^CN, HC'^N, HN'^C and DCN 
we used the collisional coefficients of the main isotopologue. We 
fixed the gas kinetic temperature and varied the density and col- 
umn density to fit the line intensities. Based on the rotation tem- 
peratures of C"^0 and CH3CN data, we assumed Tt=50 K for 
the IF and MPl positions and Tt=70 K for the MP2. Since the 
kinetic temperature is fixed in our calculations, the ratio between 
the intensities of the two lines observed is mainly dependent 
on the density while the intensity of each line is more depen- 
dent on the total column density. The obtained densities vary be- 
tween ~ a few 10^ cm"-' to ~10'' cm"-'. Our results are in agree- 
ment with Tafalla et al. (1997) who derived densities of a few 
10^ cm"^ from the CS lines and its isotopes. For some species 
we have detected only one line. In these cases we adopted a den- 
sity n(H2)=5xlO^ cm"^ in our LVG calculations for the three 
positions. 

In Table 5 we show the column densities derived towards the 
three observed positions for all the detected species. The cor- 
responding rotational diagrams are shown in appendix B (see 
Fig. Bl, B2, B3). The column densities obtained using the LVG 
approach are in agreement with the rotational diagrams calcu- 
lations which confirms our assumption that the emission is opti- 
cally thin. For those species for which the collisional coefficients 
were not known and only one transition was observed, we es- 
timated the column densities assuming optically thin emission 
and Local Thermodynamic Equilibrium with the level's popu- 
lation described by a given rotation temperature. The assumed 
rotation temperatures are based on the results of the rotational 
diagrams of other molecules with similar excitation conditions. 
For SO+, i^CN, C2D, C2S, HCO, HNCO, and C4H, we assumed 
Tro,=10 K, similar to those obtained using rotational diagrams 
for C2H and Hi^CO+. For CH3CN and CH3C2H, we adopted 
Tro,=30 K. Some cases are marked with the label "bad estimate" 
in Table 5 because the uncertainty in the column density esti- 
mates could be as high as a factor of 10. This is the case of 
CH3CN and CH3OH. The main source of uncertainty in these 
molecules is that our spectral resolution does not allow to re- 
solve all the K-components of the same rotational line. The lack 
of detection of the CH3CN 5o— >4o line towards the IF suggests 
that the emission of this transition could be self-absorbed. 

Molecular abundances were calculated assuming a canonical 
C"^0 wrt H2 abundance of 1.7x10"^. The abundance of each 
species X wrt H2 is then given by the expression fx = Nx/Nh2 = 

5. Comparison with protypical PDRs 

The 8 yum IRAC image undoubtedly shows the existence of an 
extended PDR in this region. One of the goals of this paper is to 
establish the role of the UV photons, relative to other phenom- 
ena linked to the star formation process like bipolar outflows, in 
the chemistry of this molecular cloud. For this aim, we follow a 
two-fold strategy. First of all we compare the chemical compo- 
sition in the three observed positions. Since the average kinetic 



temperature and density in the three positions are quite similar, 
differences in the chemistry can be interpreted as due to the dif- 
ferent incident UV field. Second, we compare Mon R2 with that 
in other prototypical PDRs like the Orion Bar, NGC 7023 and 
the Horsehead. 

The IF and MPl define a cut along the dense PDR in 
Mon R2, the MPl being the more shielded position. Comparing 
the molecular abundances between these two positions, we find 
that only one molecule, HCS^, has a larger (by a factor of 6) 
abundance towards the IF than towards the MPl. This suggests 
that HCS^ comes from the inner PDR like the reactive molecu- 
lar ions CO+ and HOC+ (Rizzo et al. 2003, 2005). The reactive 
ion SO^ is, however, more abundant towards the MPl. Fuente et 
al. (2003) already suggested that this molecular ion comes from 
a more shielded layer than the ions CO^ and HOC^. There are 
some species that present significantly lower abundance towards 
the IF. These species are N2H^ and HC3N proving that are eas- 
ily destroyed by UV photons and are only found in the more 
shielded gas. A peculiar chemistry is detected towards the MP2. 
This position is very rich in carbon chains such as C2H, C4H, 
C2S, and C-C3H2. As we discuss below this different chemistry 
could be due to MP2 being related with a low UV PDR. 

In Fig. 3, we compare the molecular abundances relative 
to C'**0 measured in Mon R2 with those in other prototyp- 
ical PDRs. This kind of comparison is not straightforward 
since data with different spatial resolutions and calibrations are 
used. Despite this, we can extract some qualitative conclusions. 
Mon R2 seems richer in the carbon chains HCN and CN than 
other prototypical PDRs in star forming regions such as the 
Orion Bar, NGC 7023 and the Horsehead (see Fig. 3). One could 
think that this high HCN and CN abundance could be related 
with the bipolar outflows in the region. This interpretation will 
be discussed in Sect. 5.2. The reactive ion CO^ is only detected 
in highly UV irradiated positions such as the IF in Mon R2, the 
Orion Bar and NGC 7023. The chemistry of the MP2 position in 
Mon R2 resembles that of the Horsehead, with a high C4H and 
HCO abundance. 

Abundance ratios of chemically related molecules and with 
similar excitations conditions are more accurate than absolute 
molecular abundances. On one hand, the spatial extent of the 
emission is expected the same, and the abundance ratio is not 
sensitive to the unknown beam filling factor. Moreover, as long 
as the excitation is similar for the two species, the abundance ra- 
tio is not sensitive to uncertainties in the assumed physical con- 
ditions. For these reasons, abundance ratios instead of fractional 
abundances are often used to compare with models. In Table 6 
we show some interesting abundance ratios in Mon R2 and com- 
pared them with the values in the other prototypical PDRs. 

5.1. HCN/HNC 

In this source, the HCN and HNC lines are optically thick. 
Moreover, the clear self-absorption features of these lines pre- 
vent from calculating reliable molecular abundances. Instead, we 
use the rare isotopes H^^CN and HN'^C. Using a '^C/'^C ratio of 
50 we derived X(HCN)~5xlO"^ and [HCN]/[HNC]~5-8 in the 
three observed positions. The [HCN]/[HNC] ratio is a factor ~2 
larger than in NGC 7023 and the Orion Bar (Fuente et al. 1996, 
2003). In addition to H'^CN, we have detected HC'^N towards 
the three targeted positions. The [H'^CN]/[HC'^N] is about 
10 towards the three offsets, which would imply '"*N/'^N=500 
with our assumed '^C/'-'C isotopic ratio. This ratio is similar to 
that measured by Gerin et al. (2009c) in dark clouds and low 
mass star forming regions, and consistent with the protosolar 
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Table 5. Physical parameters. 
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for the IF and MPl and Tt=70 K for the MP2 positions. When only one transition is observed, a molecular hydrogen density of 5x10' cm-'. 
Local Thermodynamic Equilibrium calculations. '' Fractional abundance derived assuming N(C'*O)/N(H2)=1.7xl0-'. '^ Assumed C'*0 wrt H2 
abundance. ^Assumed molecular hydrogen density. 
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" The density is not uniform in these regions. We indicate a representative density of the clumps/filaments. 

'' Depending on the assumed excitation conditions. 

'^ It strongly depends on the position, varying between ~1 in the shielded core to ~50 in the HCO peak. The value in the table corresponds to the 

HCO peak derived from ~14" angular resolution observations. 

'' Obtained from our '"CS data assuming '^C/'"C=50. 

References: (1) This work; (2) Rizzo et al. 2003; (3) Fuente et al. 1996; (4) Fuente et al. 2003; (5) Fuente et al. 1993; (6) Goicoechea et al. 2009; 

(7) Teyssier et al. 2004; (8) Gerin et al. 2009a; (9) Gerin et al. 2009b; (10) Goicoechea et al. 2006. 
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value of ~424. Assuming a higher '^C/'^C isotopic ratio, ~89, 
would imply '''N/'^N=890. In MPl, we have a tentative detec- 
tion of the H'^NC 7— »0 line. Using this 3-cr detection we derive 
[HC'5N]/[H15NC] = 4, similar to the values found in NGC 7023 
and the Orion Bar. A complete multitransitional study of HCN 
and its isotopologues and detailed source modelling is required 
to have a more accurate estimate of the [HCN]/[HNC] ratio and 
how it changes across the envelope in Mon R2 

5.2. CN/HCN 

Fuente et al. (1993) detected a [CN]/[HCN] abundance ratio 
larger than unity in the north-west PDR of NGC 7023 and in- 
terpreted it as the consequence of the photodissociation of HCN 
into CN. Since then, the [CN]/[HCN] ratio has been success- 
fully used as a PDR tracer in different environments, PDRs as- 
sociated with star forming regions such as the Orion Bar and 
NGC 7023, planetary nebulae, protoplanetary disks and even ex- 
ternal galaxies. In Mon R2 we use the rarer isotopes '^CN and 
H^^CN to derive the [CN]/[HCN] abundance ratio. We derive a 
[CN]/[HCN] abundance ratio ~2- 12 towards the three positions. 
This high [CN]/[HCN] abundance ratio is consistent with that 
found in other environments such as NGC 7023, the Orion Bar, 
the planetary nebula NGC 7027 and the starburst galaxy M 82. 
The [CN]/[HCN] ratio in bipolar outflow is <1 (see Bachifler & 
Perez-Gutierrez 1997). This support the interpretation of CN and 
HCN coming from the PDR instead of from the bipolar outflows. 

5.3. CS/HCS+ 

We have detected HCS^ towards the ionization front with an 
abundance of ~1.9xl0""'. This is one of the few molecules that 
shows enhanced abundance at this position and argues in fa- 
vor of HCS^ arising from the surface of the PDR. The HCS^ 
abundance measured towards the IF is a factor of~ 5 larger than 
that observed in the Horsehead nebula (Goicoechea et al. 2006) 
but comparable with that measured by Lucas & Liszt (2002) to- 
wards diffuse clouds. The HCS^ abundance towards the other 
two positions, MPl and the MP2 peaks, is lower by a factor 
of >4 and more similar to that measured in the Horsehead neb- 
ula. The derived [CS]/[HCS^] ratio is uncertain because the CS 
lines are optically thick. From our '^CS observations and as- 
suming i^cS/i^CS^SO, we obtain that the [CS]/[HCS+]~11 in 
IF and ~25 in MP2. Tafalla et al. (1997) derived column den- 
sities of N(CS)=5.6xlO'^ cm"2 towards the IF from an LVG 
analysis using the CS J=2— >7, J=5— »2 and 5—>4 lines which 
would imply [CS]/[HCS^]~7. These values are closer to the 
value [CS]/[HCS^]~10 measured in diffuse clouds than to the 
value ~175 measured towards the Horsehead nebula. 



5.4. SO^,SO, SO2, H2CS 

SO^ has been detected in warm and cold clouds with abundances 
ranging from ~10"^-10 '^ These fractional abundances can be 
explained by gas phase ion-molecule chemistry (Turner 1994). 
In PDRs, SO+ has been detected in NGC 7023 and the Orion Bar. 
The SO^ abundance in these regions is 5x10"" and 1x10"'" 
respectively. In Mon R2, we have detected SO^ in the MPl, 
with an abundance of ~5xlO"". We have not detected SO^ to- 
wards the MP2 peak. In PDRs, SO^ is primarily formed via S^ 
+ OH— >SO^ + H and removed by dissociative recombination in 
S + O. The lack of SO^ in the MP2 peak is very likely related 
to a lower abundance of OH at this position. The radical OH is 



formed by endothermic reactions in the outer layers of the PDR 
and is very sensitive to the gas kinetic temperature. 

In addition to CS, HCS+, SO and SO+, we have also detected 
the more complex sulphuretted species C2S, H2CS and SO2. SO 
has similar fractional abundance towards the three positions, but 
there are some differences in the abundances of the other com- 
pounds. The carbon chains like C2S and H2CS are more abun- 
dant towards the MP2 than towards the IF and MPl. SO2 is not 
detected towards the IF, corroborating previous interpretation of 
this molecule being easily destroyed by UV radiation (see e.g. 
Fuente et al. 2003). 

5.5. SiO 

SiO was detected in Mon R2 by Rizzo et al. (2005). They derived 
a SiO column density of 1.1 x 10" cm"^ towards the IF which is 
consistent with our upper limit of <3.0xlO" cm"^. We have de- 
tected SiO towards MPl and MP2 with an abundance of ~10"" 
which is lower than that determined by Schilke et al. (2001) to- 
wards the Orion Bar (~7xl0"" in the IF) and consistent with 
that expected in PDRs. There is no hint of enhanced SiO abun- 
dance produced by the shocks associated with the bipolar out- 
flows or the expected shock front associated to the expansion 
of the UC Hii. This is not surprising since an expansion veloc- 
ity of <1.0 km/s (Fuente et al. 2010; Pifleri et al. 2012a) is not 
high enough to destroy the core of the silicates grains. The low 
SiO abundance argues against a dominant role of shocks in the 
chemistry of this region. 

5.6. C2H, C-C3H2, C4H 

Intense C2H and C-C3H2 emission is detected in all the observed 
positions with the maximum column density towards MP2. The 
maximum fractional abundances are X(C2H)~5xlO"'' and X(c- 
C3H2)~1.7xlO""'. These values are a factor of ~3 lower than 
those measured towards the Horsehead and the Orion Bar but 
similar to those measured in NGC 7023. Taking into account 
the uncertainties in the calculation of the fractional abundances, 
these differences are not significant. Note that for the Horsehead 
we have taken the single-dish column densities published by 
Teyssier et al. (2004) to compare with our data because interfero- 
metric observations filter out the extended emission tracing only 
the small spatial scales. The [c-C3H2]/[C2H] ratio is quite uni- 
form with a value of ~0.03. The large chain C4H had only been 
detected in low-UV PDRs (IC 63, p-Oph and the Horsehead) 
thus far (Teyssier et al. 2004; Pety et al. 2005). The [C4H]/[C2H] 
ratio is maximum at the MP2 confirming that the PDR at this po- 
sition is particularly rich in carbon chains. 

5.7. HCO 

The formyl radical, HCO, was detected in the interstel- 
lar medium by Snyder et al. (1976). Pioneering studies by 
Schenewerk et al. (1988) pointed out to an association of large 
HCO abundances with the gas around Hii regions. Schilke et al. 
(2001) detected HCO towards NGC 2023, the offset (-30",80") 
in NGC 7023, the position referred to as Bar(CO) in the Orion 
Bar and the IF in S140. From a limited mapping they con- 
cluded that HCO is associated with the PDR component and ex- 
plained the enhancement of HCO abundance as the consequence 
of the photodissociation of the H2CO molecules released from 
the grain mantles by UV photons. Garcia-Burillo et al. (2002) 
mapped this radical in the starburst M 82 and also interpreted it 
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Fig. 3. Comparison of the abundance of the most representative species among the prototypical PDRs: Mon R2, NGC 7023, Orion 
Bar and the Horsehead. 



as arising in the giant PDR in the disk of this galaxy. Gerin et 
al. (2009a) carried out a high angular resolution mapping of the 
Horsehead in this radical. The maximum of the HCO emission 
was found towards a position where the PAH band emission is 
particularly intense. The HCO emission becomes fainter in the 
more shielded molecular gas. After updating some reaction rates, 
in particular O+CH2 — > HCO + H, they found that gas-phase 
chemistry could account for the high HCO abundance without 
invoking the photodesorption of H2CO from the ice mantles. 
In Mon R2, we have detected HCO towards MP2 with a frac- 
tional abundance of ~1.7xl0 '". Towards the IF and the MPl, 



the HCO abundance is at least a factor of ~2 lower than towards 
the MP2. Note that Schilke et al. (2001) did not detect HCO to- 
wards the IF but only in the position referred to as Bar(CO) in 
their paper that could correspond to a lower UV radiation field. 



5.8. H2CO 

We have derived a H2CO abundance of a;l-2xl0 ** in Mon R2. 
This abundance is close to that found in hot cores (see e.g. Sutton 
et al 1995) and higher than that expected in a PDR. Guzman et 
al. (201 1) studied the formation of H2CO in the Horsehead. They 
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derived an H2CO abundance of 4-6x10'*^ in the PDR peak, 
similar to that found in the shielded core position. They could not 
explain the H2CO abundance in this PDR only with gas phase 
chemistry. However, they obtained a good agreement between 
model predictions and observations with a gas phase+surface 
chemistry model in which H2CO can be formed on the grain 
surfaces and photo-desorbed to the gas phase, in addition to the 
gas phase route. In Mon R2, the emission of H2CO is ubiquitous 
(Giannakopoulou et al. 1997) and the low spatial resolution of 
our observations makes it difficult to discern its origin: the PDR, 
the more shielded molecular gas or the bipolar outflows. 

5.9. Complex molecules: CH3OH, CH^CN, HC3N 

In addition to PDR-like molecules, we have detected in our sur- 
vey other complex species that are not expected to be abun- 
dant in PDRs and are usually associated with the chemistry of 
warm clouds. This is the case of CH3OH, CH3CN and HC3N. 
In the case of methanol, our fractional abundance estimates are 
very uncertain. The abundance of CH3CN is also quite constant 
and equal to ^7x10 ". For HC3N we measured abundances of 
7xl0"^'-10"'". These values are lower than those found in hot 
cores, but very similar to those found in warm molecular clouds 
such as the Orion Ridge and in the envelopes of massive proto- 
stars. 

Complex molecules are rapidly destroyed by UV radiation. 
Assuming a typical photodissociation rate of 10"'° s"', the life- 
time of complex molecules would be >10^ yr only for Av>10 
mag. Fuente et al. (2010) and Pilleri et al. (2012a) estimated an 
expansion velocity of < 1 .0 km s ' for this UC Hii which implies 
a dynamical time of ~10^ yr. This means that complex molecules 
of the initial core could have survived as long as they were within 
clumps in which the visual extinction is >10 mag. High spa- 
tial resolution observations are required to determine the region 
from which the emission of complex molecules arises and put 
additional constraints to their formation mechanism. 



5.10. Deuteration: DCN, C2D 

We have detected DCN and C2D towards the MPl and 
MP2. The DCN abundances are X(DCN)=1.0xlO-''' (MPl) 
and X(DCN)=4.2xlO-" (MP2). The C2D abundances are 
X(C2D)=2.3xlO-"' (MPl) and X(C2D)=2.6xlO-i" (MP2). 
These abundances correspond to deuterium fractionation values 
of -0.03 -0.05 for both HCN and C2H. The deuterium fraction- 
ation of HCN is a factor of ~10 larger than that measured by 
Parise et al. (2009) in the position Bar(HCN) of the Orion Bar. 
The deuterated species C2D was not detected in Orion with a 
lower limit to the C2D/C2H ratio of <0.04. Pety et al. (2007) 
detected DCO^ in a cold clump (Tj^~ 10-20 K) towards the 
Horsehead and measured [DCO+]/[H'^CO+] =0.02. There is not 
any intense transition of DCO^ in the frequency range covered 
by our survey. Therefore, we could not obtain an estimate of the 
DCO^ abundance. 

Several mechanisms have been proposed to explain the high 
values of deuterium fractionation (higher than the deuterium 
abundance in the Universe, D/H~10"^) observed in the in- 
terstellar medium. Deuterated isotopologues of methanol and 
formaldehyde have been detected in hot cores and corinos (see 
e.g. Parise et al. 2002, 2004; Fuente et al. 2005) where the gas 
kinetic temperatures is >100 K. In these cases, the fractionation 
is thought to occur on grain surfaces. The deuterated compounds 
are released to the gas phase when the ice is evaporated, produc- 



ing a transient deuterium enhancement in gas phase. In Mon R2, 
we have detected the deuterated compounds of HCN and C2H. 
Since these molecules are not mainly formed on grain surfaces, 
it is more plausible to interpret that the deuteration has occurred 
in gas phase. 

In molecular clouds, deuterium is mainly locked into HD. 
Efficient transfer of deuterium from this reservoir to other 
species occurs by means of ion-molecule reactions. In cold 
clumps, Tj^»10-20 K, deuteration is usually transferred via re- 
actions with H2D^. This mechanism is very efficient, for in- 
stance, in pre-stellar cores where molecules are highly depleted 
(see e.g. Caselli et al. 2003). For slightly higher temperatures of 
Tj^»30-50 K, the transfer is more efficient via CH2D^ (Roueff 
et al. 2007, Parise et al. 2009). The high temperatures measured 
in Mon R2 favor a deuteration mechanism based on CH2D^. 
Our measurements are in agreement with gas-phase model pre- 
dictions by Roueff et al. (2007) and Parise et al. (2009) for 
Tj^~30-50 K and densities ~a few 10^ cm"^. These densities are 
not unrealistic in this region (Rizzo et al. 2005). One possibility 
is that C2D comes from the densest part of this PDR. Another 
possibility is that a fraction of these molecules are formed on 
grain surfaces or their chemistry is related with some of the evap- 
orated species. 

6. Chemical diagnostics 

From the discussion in Sect. 5 it seems clear that the emission 
of some of the detected species comes from the extended PDR 
around this Hii region. In the following we try to understand 
the differences between MPl and MP2 in terms of PDR chem- 
istry. We have run a grid of isochoric models using the updated 
version (1.4.2) of the PDR Meudon code (Le Petit et al. 2006, 
Goicoechea et al. 2007) and the parameters listed in Table 7. Our 
grid of models range in density between 10"* cm"^ and 10^ cm"^, 
and in FUV field between 10 and lO*" Habing fields. We have 
run the same grid of models using the standard galactic and the 
Orion extinction curves but the results are essentially the same 
for the molecular abundance ratios discussed below. For each 
model we represent the cumulative column densities summed 
up to 10 mag. For higher visual extinctions, UV radiation has 
a negligible effect on the physical conditions and chemistry of 
the molecular gas. Our goal is to investigate the influence of the 
incident UV radiation and the density on the studied molecu- 
lar abundance ratios and explore the possibility of using them 
as chemical diagnostics. These models are useful to interpret 
molecular observations, although to reproduce the results on a 
particular PDR it is required a good knowledge of its geometry. 
Some representative cumulative column density ratios have 
been selected as candidates for chemical diagnostics. Fig. 4 
shows these ratios as a function of the incident UV field and the 
hydrogen nuclei density for the standard interstellar extinction 
curve. Since we are using a gas phase model, we have avoided 
those species for which the formation on grain surfaces could 
be important. This is the case of H2CO, CH3CN and CH3OH. 
Another important reason to neglect these complex molecules is 
that they could come from a shielded component of dense gas 
instead of the PDRs. We are aware, however, that the injection 
of complex molecules in gas phase by photo-desorption changes 
the gas chemical composition and could influence the chemistry 
of more simple species. Consider, for instance, HCO, which is 
the photodissociation product of H2CO. The abundance of HCO 
in the PDR would increase if the evaporation of H2CO from 
the grain surfaces were included improving the agreement with 
the observations. As discussed below, PAHs and grain destruc- 
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Table 7. Input parameters for the Meudon PDR code. 



Parameter 






Value 


Go 


Radiation field intensity (Habing) 


10 to 10" 


G«' 


External Radiation field 
(Habing) 


intensity 


1 


n 


Hydrogen nuclei density 




Iff* to 10' cm-' 


Ay 


Cloud depth 




10 mag 


Extinction 


Standard Galactic 






Rv 


Av/E{B - V) 




3.1 


i 


Cosmic ray ionization rate 




5x10-" S-' 


^min 


dust minimum radius 




3x10-' cm 


^max 


dust maximum radius 




3x10-^ cm 


a 


MRN dust size distribution index 


3.5 


He/H 


Helium abundance 




0.1 


O/H 


Oxygen abundance 




3.2x10-" 


C/H 


Carbon abundance 




1.3x10-* 


N/H 


Nitrogen abundance 




7.5x10-5 


S/H 


Sulphur abundance 




1.9x10-5 


Fe/H 


Metal abundance 




1.5x10-** 



lion are also required to account for the observed abundances of 
small hydrocarbons. 

Our chemical modeling shows that for the range of param- 
eters considered, the CN/HCN ratio is not dependent on the 
incident UV field but on the density. Therefore, the CN/HCN 
ratio is a good probe of the existence of PDRs but it is not 
a good diagnostic to differentiate between PDRs with differ- 
ent UV fields for Go>100. The value measured in Mon R2, 
CN/HCN=2-12, is consistent with a hydrogen nuclei density 
of »10^ cm-^ (n(H2)~5xlO^ cm"^) in agreement with our den- 
sity estimates. As commented in Sect. 5, the HCN and CN 
abundances in this region are larger than in other prototypical 
PDRs. We cannot discard that the existence of bipolar outflows 
could contribute to enhance the abundance of these molecules. 
However, the measured CN/HCN=2-12 is more consistent with 
a PDR origin. 

According with our models, the HCO^/HCN ratio must be 
a good tracer of the UV field and density. Fig. 4b shows that 
the HCO^/HCN ratio is lower than 1 for a wide range of phys- 
ical conditions but increases to larger values for high values of 
Go (>10-') and of the density (~10*' cm--'). Taking into account 
the extreme values of Go and n in Mon R2, one would expect a 
higher ratio than that observed. As mentioned above, the detailed 
geometry of the PDR has to be taken into account when compar- 
ing with chemical models. In Mon R2 several gas components 
lie along the line of sight. The dense PDR is surrounded by a 
lower density envelope that very likely contributes to the emis- 
sion of the low J rotational lines of HCO^ and HCN. Although 
we use the rare isotopologues, H'^CO^ and H'^CN, to compute 
this abundance ratio, optical depth effects could still be impor- 
tant. 

One of the best diagnostic of the UV incident field is the 
CO+/HCO+ ratio. This ratio is >0.05 only when flie UV field 
is higher than 10^. This is consistent with the lack of detec- 
tion of CO^ in the Horsehead nebula and in the MP2 position of 
Mon R2. The same dependency on the UV incident field stands 
also for the HOC^/HCO^ ratio. On the contrary, the formyl rad- 
ical is abundant only in PDRs with moderate or low UV radia- 
tion field (see also Gerin et al. 2009ab). The HCO/HCO+ ratio 
takes values >0.2 only for UV fields <10l The CO+/HCO+ and 
HCO/HCO^ abundance ratios are excellent diagnostics to differ- 
entiate between low and high UV PDRs. There is a good qualita- 



tive agreement between models and observations although they 
still fail to quantitatively predict these ratios. 

Observations show that the SO^/SO ratio is also higher for 
dense and highly ionized regions and the value towards Mon R2 
is well predicted by our gas-phase chemical models. However, 
observationally this ratio does not follow the same trend that the 
CO^/HCO^ (see e.g. Fuente et al. 2003). This ratio is influenced 
by the amount of sulphur in gas phase that could be changed, for 
instance, by the existence of bipolar outflows and slow shocks. 
Therefore, we do not propose it as a good diagnostic of PDRs. 

Small hydrocarbons can be good tracers of PDRs, and they 
are specially abundant in low UV PDRs. According with our gas 
phase calculations, the C-C3H2/C2H ratio presents a fairly uni- 
form value in a wide region of the parameter space. Although 
our models can explain the value measured in MP2, they fall 
short to account for the value measured in the Horsehead. In the 
case of the C4H/C2H ratio, our models fall short to reproduce 
the observed values in both, MP2 and the Horsehead. Different 
authors have pointed out that the observed abundances of these 
hydrocarbons in PDRs requires of another production mecha- 
nism than gas-phase reactions, possibly the destruction of PAHs 
and/or very small grains (VSGs). It is a matter of a lower value 
of the UV field that prevents too fast destruction of hydrocar- 
bons, but having destruction of PAHs/VSGs at a reasonable rate. 
Going further requires a full modelling of formation/destruction 
processes and it is beyond the scope of this paper. 

Summarizing, the CO+/HCO+ and HCO/HCO+ ratios are 
excellent chemical diagnostic of the UV field, while the 
CN/HCN and HCN/HNC ratios are more dependent on the den- 
sity. Applying these chemical diagnostics to the three positions 
observed in Mon R2, we can conclude that the incident UV ra- 
diation field in the IF and MPl are different, being the highest 
in IF, with a value consistent with that estimated by Rizzo et al. 
(2003). Towards the MP2, the incident UV field must be lower. 
The morphology of the 8 //m emission (see Fig. 1) suggests that 
this PDR is the consequence of the border of the molecular cloud 
being illuminated by IRS 1 . 

7. Summary and Conclusions 

We have carried out a mm survey towards three positions of 
the PDR that represent different physical and chemical environ- 
ments: (i) the ionization front (IF); two positions in the molec- 
ular cloud, (ii) MPl and (iii) MP2. Our goal has been to inves- 
tigate the chemistry of the molecular gas around the UC Hii re- 
gion and the possible variations due to the different local physi- 
cal conditions. Our results can be summarized as follows: 

- We have detected more than thirty different species (includ- 
ing isotopologues and deuterated compounds) in Mon R2. 
In particular, we have detected SO^ and C4H towards this 
region, which are well known tracers of PDRs. In addition 
to SO^ and C4H, the list of identified species includes typi- 
cal tracers of PDRs like CN, HCN,HCO, C2H, and C-C3H2 
but also other complex molecules more common in warm 
molecular clouds such as CH3CN, H2CO, HC3N, CH3OH 
or CH3C2H. The origin of these complex species, dense and 
well shielded clumps within the molecular cloud or PDR, is 
still to be investigated. The existence of high velocity molec- 
ular outflows in the region could also affect the abundance of 
some molecules such as HCN and the sulfuretted species. 

- Within the PDR component, the comparison of the fractional 
abundances measured in Mon R2 with those in the prototyp- 
ical PDRs shows that the positions IF and MPl present a 
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Fig. 4. Cumulative column density ratios in a plane slab of Av=10 mag illuminated by the left side for a grid of UV fields and 
hydrogen nuclei densities. The ratios have been calculated using the Meudon code (Le Petit et al. 2006, Goicoechea et al. 2007). 
Go is given in unit of the Habing field. Levels drawn in red correspond to the values measured in Mon R2. Symbols indicate the 
physical conditions of the protypical PDRs: Mon R2 (square), Orion Bar (triangle), NGC 7023 (circle) and Horsehead (cross). Note 
that the CO+/HCO+ and HCO/HCO+ ratios are excellent diagnostics of the UV incident field. The IF and MP positions in Mon R2 
present abundances characteristics of dense PDRs with high UV field while the chemistry of the MP2 position is more typical of a 
low UV PDR. 



chemistry similar to that in high-UV PDRs (Go>10^ Habing 
fields) while the chemistry in the position MP2 resembles 
that of the Horsehead ( <10^ Habing field). Chemical mod- 
els predict that [CO+]/[HCO+] and [HCO]/[HCO+] ratios are 
good diagnostics to differentiate between these two types of 
PDRs. 
- The deuterated species DCN and C2D are detected in our 
spectral survey. The observed deuteration fractionations, 
[DCN]/[HCN]~0.03 and [C2D]/[C2H]~0.05, are among the 
highest in warm regions. 

One important question is whether these stationary models 
are adequate for PDRs associated with rapidly evolving star 
forming regions. The time to reach the equilibrium for each 
species at a given Ay depends on the density and the UV ra- 
diation (see e.g Bayet et al. 2009). The higher the density and 



incident UV radiation, the fastest the equilibrium is reached at a 
given A,,. For the physical conditions in Mon R2, simple species 
are expected to have reached the equilibrium in ~10^ yr which 
explains the success of our chemical diagnostics. However, we 
cannot exclude small deviations from equilibrium especially in 
those regions at high extinction and with low incident UV field. 
In following works, we will investigate the role of grain surface 
chemistry and time dependent effects on the chemistry of the 
molecular gas in Mon R2. 
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Appendix A: Identified lines 



Table A.l. Observed transitions and line intensities 



Frequency 


Transition 


Intensity (IF) 


Intensity (MPl) 


Intensity (MP2) 


(MHz) 




(Kxkms"') 


(Kxkms-') 


(Kxkms-') 


83802 


H(66)<5 


1.23(0.11) 


0.24(0.06) 




84107 


C(40)a'+He(40)Q' in the image band 




0.66(0.14) 




84521 


CH30H5_i,5-^4o,4 


0.48(0.08) 


0.35(0.04) 


0.62(0.04) 


84914 


H(60)7 


2.26(0.12) 


0.70(0.04) 




85197 


CS 2— >i in the image band 




0.39(0.19) 


0.39(0.09) 


85339 


C-C3H2 2i,2^io,l 


1.18(0.20)" 


1.50(0.07)° 


2.09(0.09)" 


85348 


HCS+ 2^1 


0.76(0.09) 


0.32(0.02) 


0.33(0.09) 


85457 


CH3C2HJ^^4tk=0,l 


0.41(0.07) 


0.47(0.16) 


0.53(0.07) 


85634 


C4H9^SJ=79/2^i7/2 






0.29(0.08) 


85673 


C^n 9^8^=17/2^15/2 






0.53(0.15) 


85688 


H(42)a 


15.97(0.13) 


4.08(0.13) 


0.89(0.16) 


85731 


C(42)a 




0.23(0.06) 




85925 


NHzDiu-io.i' 




0.19(0.04) 




86055 


HC'^N 1-^0 


0.29(0.07) 


0.52(0.14) 


0.37(0.06) 


86094 


SO 22^7, 


0.38(0.05) 


0.31(0.12) 


0.54(0.01) 


86339 


H"CN7^0 


1.78(0.08) 


2.11(0.03) 


1.64(0.03) 


86488 


H(70)6 


0.95(0.12) 






86671 


HCO loi-^Ooo }=3/2^1/2 F=2^1 




0.22(0.05) 


0.54(0.05) 


86690 


H(74)^ 


0.99(0.17) 






86708 


HCO io.i^Oo.o }=3/2^1/2 F=1^0 




0.15(0.03) 


0.37(0.05) 


86754 


H"CO+ 1^0 


0.55(0.09) 


1.43(0.02) 


1.14(0.09) 


86777 


HCO ioi^Ooo ]=l/2^1/2 F=l^l 




0.12(0.03) 


0.31(0.04) 


86847 


SiO 2^1 




0.32(0.04) 


0.35(0.01) 


87091 


HN"C i^O 


0.31(0.06) 


0.57(0.04) 


0.35(0.04) 


87284 


C2H N=7^0 ]=3/2^1/2 F=i^i 


0.42(0.08) 


0.68(0.02) 


0.69(0.04) 


87317 


C2H N=1^0 J=3/2^1/2 7=2^1 


4.37(0.38) 


5.47(0.49) 


5.77(0.46) 


87329 


C2H N=7^0 i=3/2^1/2 F=1^0 


2.12(0.11) 


2.97(0.06) 


3.00(0.36) 


87402 


C2H N=7^0 ]=l/2^1/2 F=l^l 


2.17(0.14) 


3.24(0.06) 


3.38(0.13) 


87407 


C2H N=i^O J=l/2^1/2 F=0-^1 


0.84(0.15) 


1.15(0.05) 


1.10(0.06) 


87446 


C2H N=i^O ^=1/2^1/2 F=1^0 




0.77(0.06) 


0.68(0.05) 


87615 


H(65)(5 


1.86(0.08) 


0.41(0.09) 




87706 


Unidentified (F,„,„j,= 110433 MHz) 


1.47(0.30) 






87938 


"CO 1^0 in the image band 


0.96(0.10) 


1.05(0.05) 


0.84(0.05) 


88358 


C'**0 7->0 in the image band 




0.36(0.03) 


0.27(0.04) 


88406 


H(52)/3 


6.77(0.08) 


1.01(0.07) 


0.47(0.06) 


88578 


Unidentified (F™<,g,= 109560 MHz) 


0.47(0.06) 






88632 


HCN 1-^0 F=l-1, F=2-l and F=0-1 


8.47(0.06) 


20.58(0.05) 


23.24(0.04) 


88866 


H'^NC i^O' 




0.14(0.04) 




88942 


Unidentified (F™<,g,=109196 MHz) 






0.20(0.04) 


89112 


Unidentified (F™„j,= 109027 MHz) 


0.24(0.05) 






89165 


Unidentified (F™<,g,= 108974.5 MHz) 


0.66(0.02) 






89189 


HCO+ 1^0 


11.59(0.02)" 


13.90(0.22) 


17.67(0.11)" 


89199 


H(59)r 


2.50(0.21) 






90174 


H(69)e and H(73)f 


1.63(0.12) 






90664 


HNC 1^0 


3.19(0.05) 


10.58(0.14) 


8.33(0.06) 


90979 


HC3N 10^9 


0.90(0.06) 


2.41(0.08) 


1.41(0.03) 


91406 


H(39)a in the image band 


0.61(0.13) 






91663 


H(64)^ 


1.86(0.11) 


0.52(0.08) 




91942 


CH3'3CN5o^4o 


0.59(0.16) 






91971 


CH3CN 53^43 


0.64(0.14) 


0.11(0.05) 


0.11(0.03) 


91980 


CH3CN 52^42 






0.12(0.03) 


91987 


CH3CN5^^4^k=0,l 




0.52(0.05) 


0.41(0.04) 



' The profile shows self-absorption; ' Tentative detection; obtained by degrading the HERA maps to the angular resolution of the 3mm lines. 
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Table A.l. Observed transitions and line intensities (continuation) 



Frequency 


Transition 


Intensity (IF) 


Intensity (MPl) 


Intensity (MP2) 


(MHz) 




(Kxkms"') 


(Kxkms"') 


(Kxkms"') 


92034 


H(41)a 


20.45(0.17) 


4.01(0.18) 


1.39(0.05) 


92080 


C(41)a 


0.38(0.08)' 






92436 


Unidentified (F,„,„j,= 105702 MHz) 


0.39(0.05) 






92494 


"CS 2-»i 


0.30(0.05) 


0.42(0.06) 


0.37(0.03) 


93174 


N2H+ 1^0 ¥=1^1 , 1^2 and 1^0 


1.53(0.06) 


8.66(0.09) 


2.47(0.04) 


93607 


H(51)/3 


6.04(0.10) 


1.06(0.08) 


0.35(0.08) 


93776 


H(58)r 


2.77(0.17) 


0.56(0.09) 




93827 


H(72)^ 


0.89(0.09) 






93870 


C2S 78^6, 


0.14(0.04) 




0.16(0.02) 


94073 


H(68)6 


1.03(0.07) 






103252 


Unidentified (F™„g,=94889 MHz) 


0.52(0.04) 






103266 


C4Hv=l 21/2-^19/2' 


0.47(0.05) 






103915 


H(56)r 


3.02(0.07) 


0.46(0.07) 




104029 


SO2 5i.3^2o,2' 




0.16(0.06) 


0.11(0.03) 


104062 


Unidentifed (F,™^,=93900 MHz) 






0.16(0.03) 


104239 


SO2 iOi,,^70o,io' 






0.10(0.03) 


104533 


H(51)/3 in tire image band 


0.44(0.03) 






104617 


rl2Co J[ 2 — *2|j 


0.19(0.03) 


0.38(0.03) 


0.21(0.02) 


104667 


C4iill^l0]=23/2^21/2 




0.19(0.03) 


0.30(0.06) 


104705 


C4H 11-^10 ]=21/2^19/2 




0.27(0.05) 


0.30(0.04) 


104966 


N2H^ l—>0 in the image band 




0.52(0.03) 




105302 


H(49)/3 


6.22(0.10) 


0.66(0.06) 


0.39(0.05) 


105410 


H(61)(5 


1.70(0.07) 






105516 


Unidentified (Fi,„„j,= 92625 MHz) 


0.36(0.10) 






106079 


H(69)^ 


0.71(0.09) 






106104 


H(41)a in the image band 


0.45(0.10) 






106348 


C2S N=8^7 J=9^8 






0.13(0.02) 


106737 


H(39)a 


18.80(0.11) 


2.33(0.07) 


0.57(0.06) 


107014 


CH30H5i,3^4o,4 




0.14(0.06) 


0.55(0.07) 


107206 


H(65)e 


1.23(0.12) 






107254 


Unidentified (Fi,„<,g,=90887 MHz) 


0.59(0.18) 






107424 


Unidentified (Fi,„„j,=90720 MHz) 


0.51(0.03) 






107476 


HNC i— >0 in the image band 


0.25(0.07) 


0.46(0.03) 


0.48(0.04) 


108427 


"CNN=7^0J=7/2^7/2 




0.27(0.01) 


0.16(0.01) 


108780 


"CN N=7^0 J=i/2^7/2 




0.24(0.01) 


0.39(0.03) 


108894 


CH3OH Oo,o^i-i.i 




0.20(0.09) 


0.19(0.03) 


108950 


HCO+ i->0 in the image band 


0.21(0.21) 


0.47(0.13) 


0.72(0.05) 


109174 


HC3N 12^11 


0.70(0.04) 


2.41(0.04) 


1.45(0.04) 


109206 


Unidentified (F™„j,=88935 MHz) 






0.13(0.03) 


109252 


S0 23^i2 


0.21(0.18) 


0.34(0.01) 


0.59(0.01) 


109536 


H(55)r 


3.25(0.10) 






109782 


C'^O 1^0 


3.47(0.06) 


4.63(0.05) 


3.75(0.06) 


109906 


HNCO 5^4 




0.13(0.04) 




109998 


Unidentified (Fi,„„j,=88142 MHz) 


0.43(0.05) 






110025 


C'^N N=1^0J=3/2^1/2 F=2^1 and 1^0 




0.13(0.03) 


0.25(0.04) 


110119 


Unidentified (F™<,g,=88022 MHz) 






0.18(0.04) 


110201 


'^CO 1^0 


42.31(0.03) 


49.18(1.07) 


41.61(0.15) 


110364 


CH3CN 63^53 






0.21(0.06) 


110375 


CH3CN 62^52 




0.12(0.04) 


0.13(0.04) 


110383 


CHsCN 6k^5k k=0,l 


1.13(0.13) 


0.59(0.04) 


0.34(0.05) 


110601 


H(60)(5 


1.66(0.13) 






205081 


CH3C2H 12o^11q 


0.49(0.08) 


0.71(0.10) 




205292 


C'*0 2^1 in the image band 






1.01(0.08) 


205986 


H2CS 6o,6^5o,5 




0.26(0.06) 





' The profile shows self-absorption; ' Tentative detection; obtained by degrading the HERA maps to the angular resolution of the 3mm lines. 
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Table A.l. Observed transitions and line intensities (continuation) 



Frequency 


Transition 


Intensity (IF) 


Intensity (MPl) 


Intensity (MP2) 


(MHz) 




(Kxkms"') 


(Kxkms"') 


(Kxkms"') 


206176 


SO 45^34 


1.41(0.09) 


1.30(0.07) 


1.60(0.09) 


206676 


Unidentified (F™<,g,=218175 MHz) 


1.31(0.12) 


0.79(0.06) 


0.65(0.05) 


206907 


Unidentified (Fi™j,=217944 MHz) 


0.49(0.11) 






207836 


Unidentified (F™<,j,=217015 MHz) 






0.40(0.08) 


208700 


SO2 52,2^2i,, 






0.35(0.07) 


208767 


Unidentified (F™<,g,=221286 MHz) 




0.29(0.07) 




209097 


Unidentified F,„„j,=220956 MHz) 


0.64(0.06) 


0.76(0.05) 


0.86(0.12) 


209201 


H2CS6i5-^5i4 


0.39(0.04) 


0.53(0.08) 


0.70(0.16) 


209230 


HC3N 23-^22 


0.22(0.10) 


0.28(0.11) 


0.59(0.15) 


209277 


Unidentified (Fi,„„j,=220777 MHz) 


0.37(0.08) 


0.32(0.08) 




209591 


Unidentified (F™<,g,= 220462 MHz) 


1.84(0.54) 






209655 


'^CO 2— >i in the image band 


1.97(0.10) 


2.13(0.07) 


2.18(0.11) 


209894 


H(44)7 


2.37(0.16) 






209999 


C(44)7 






0.51(0.11) 


210495 


C'*0 2->i in the image band 


1.18(0.07) 


0.32(0.06) 


1.28(0.09) 


210502 


H(31)a 


15.18(0.20) 






216373 


CCD N=i^2 J= 7/2^5/2 




0.31(0.05) 


0.49(0.10) 


216428 


CCD N=i^2 J=J/2^i/2 




0.31(0.05) 


0.40(0.07) 


216711 


H2S 22,0— >2i_i 






0.21(0.05) 


217105 


SiO 5^4 




0.22(0.05) 


0.23(0.01) 


217239 


DCN J=5->2 


1.23(0.11) 


2.94(0.06) 


1.32(0.06) 


217467 


'^CN N=2^i J=J/2^J/2 ¥=4^3, 3^2 and 2^1 




0.12(0.01) 


0.38(0.03) 


217822 


C-C3H2 6i_6^5o,5 and 6o.6-»-5i,5 


0.80(0.20) 


1.20(0.20) 


1.86(0.20) 


217835 


CO 2— >i in the image band 


1.80(0.16) 


1.95(0.45) 


3.04(0.40) 


217940 


C-C3H2 5,,4^42,3 




0.52(0.05) 


1.13(0.05) 


218160 


C-C3H2 52.4^4i.3 




0.33(0.15) 


0.40(0.03) 


218182 


Unidentified (F,„,„j,=230192 MHz) 




0.50(0.09) 


0.58(0.11) 


218222 


H2COJo^-^2o2 


3.78(0.13) 


7.71(0.10) 


6.24(0.12) 


218325 


HC3N 24^23 


0.17(0.06) 


0.27(0.04) 


0.32(0.03) 


218440 


CH30H42 2^-?12 


0.63(0.13) 


0.74(0.04) 


0.35(0.06) 


218476 


H2CO 52.2^22.1 


1.02(0.12) 


1.77(0.09) 


1.42(0.07) 


218667 


Unidentified (F,„„j,=229707 MHz) 


0.90(0.11) 






218760 


H2CO 52,1-^22.0 


0.87(0.08) 


1.75(0.05) 


1.41(0.06) 


219560 


C'*0 2^i 


10.68(0.17) 


14.0(0.07) 


12.74(0.07) 


219597 


Unidentified (Fi,„„j,=228776 MHz) 


1.90(0.13) 






219852 


H2CN5o,3^2o,2' 


1.03(0.29) 


0.57(0.05) 


0.48(0.06) 


219911 


H2"C0 5i2-2ii' 


1.06(0.15) 






219949 


SO 65-»J4 


1.45(0.07) 


2.07(0.01) 


1.83(0-03) 


259012 


H"CN 3^2 * 


1.65(0.21) 


2.50(0.25) 


1.10(0.06) 


260255 


H"CO+ 3^2 * 


2.00(0.04) 


2.12(0.04) 


1.52(0.03) 


262004 


C2HN=5-^2J= 7/2^5/2* 


14.59(0.60) 


15.60(0.60) 


12.89(0.47) 



" The profile shows self-absorption; ' Tentative detection; obtained by degrading the HERA maps to the angular resolution of the 3mm lines. 
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Table A.2. Summary of recombination lines 







IF (0" 


,0") 






Freq(MHz) 


Line 


I(Klflns-') 


v,.,r(kms ') 


Av(kms-') 


Tmb(K) 


83801.829 


H(66)(5 


1.23(0.11) 


13.3(1.3) 


27.8 (2.9) 


0.04 


84914.394 


H(60)7 


2.26(0.12) 


11.8(0.8) 


30.7(1.8) 


0.07 


85688.390 


H(42)a 


15.97(0.13) 


12.4(0.1) 


30.9 (0.28) 


0.48 


85731.144 


C(42)a 


0.23 ( 0.06) 


12.2(1.6) 


9.0 ( 2.9) 


0.02 


86488.417 


H(70)e 


0.95 (0.12) 


14.8 (2.6) 


36.9 (5.5) 


0.02 


86689.621 


H(74)^ 


0.99(0.17) 


8.2 (2.7) 


32.5 (8.0) 


0.03 


87614.996 


H(65)(5 


1.86(0.08) 


9.7 (0.7) 


30.3 (1.5) 


0.06 


88405.687 


H(52)/3 


6.77 (0.08) 


9.3 (0.2) 


34.2 (0.5) 


0.19 


89198.545 


H(59)r 


2.50(0.21) 


7.2(1.1)" 


26.8 (2.5) 


0.09 


90164.065 


H(73K* 


0.63(0.10) 


9.8(2.3) 


31.0(7.0) 


0.02 


90174.339 


H(69)6* 


0.98(0.12) 


8.4(1.3) 


32.0 (3.0) 


0.03 


91663.133 


H(64)(5 


1.86(0.11) 


9.0(1.0) 


33.2 (2.0) 


0.05 


92034.434 


H(41)q' 


20.45 (0.17) 


9.6 (0.2) 


30.8 (0.3) 


0.62 


92080.355 


C(41)a' 


0.38 (0.08) 


7.0(1.2) 


11.7(3.1) 


0.03 


93607.316 


H(51)/3 


6.04(0.10) 


9.8 (0.3) 


30.9 (0.6) 


0.18 


93775.871 


H(58)r 


2.77(0.17) 


8.7 (0.9) 


30.5 (2.1) 


0.08 


93826.790 


H(72K 


0.89 (0.09) 


10.7 (1.5) 


29.6 (3.8) 


0.03 


94072.811 


H(68)6 


1.03 (0.07) 


11.4(0.9) 


26.8(2.1) 


0.04 


103914.838 


H(56)r 


3.02 (0.07) 


12.6 (0.4) 


33.0 (0.9) 


0.09 


105301.857 


H(49)/3 


6.22(0.10) 


12.4 (0.3) 


31.5 (0.6) 


0.19 


105410.216 


H(61)<5 


1.70(0.07) 


12.7 (0.7) 


29.5 (1.4) 


0.05 


106079.529 


H(69K 


0.71 (0.09)) 


14.0(1.9) 


28.8(4.1) 


0.02 


106737.357 


H(39)q' 


18.80(0.11) 


12.53 (0.09) 


31.0(0.2) 


0.57 


107206.108 


H(65)e 


1.23(0.12) 


10.31 (1.8) 


35.5 (3.9) 


0.03 


109536.001 


H(55)7 


3.25(0.10) 


12.5 (0.5) 


33.5 (1.2) 


0.09 


110600.675 


H(60)(5 


1.66(0.13) 


10.4(1.3) 


34.0 (3.2) 


0.04 


209894.050 


H(44)7 


2.37(0.16) 


13.3(1.1) 


31.6(2.6) 


0.07 


210501.771 


H(31)q'^ 


15.18 (0.20) 


13.7 (0.3) 


33.0 (0.4) 


0.43 






MPl (15' 


',-15") 






Freq(MHz) 


Line 


I(K]flns-') 


v,„(kms ') 


Av(kms-') 


Tmb(K) 


83801.829 


H(66)<5 


0.24 (0.06) 


16.1(2.4) 


17.0 (4.9) 


0.01 


84914.394 


H(60)r 


0.70 (0.04) 


15.8 (0.9) 


31.5(1.8) 


0.02 


85688.390 


H(42)a 


4.08(0.13) 


15.3 (0.5) 


29.7(1.2) 


0.13 


85731.144 


C(42)a 


0.22 (0.03) 


8.7 (0.7) 


12.0(1.7) 


0.02 


87614.996 


H(65)(5 


0.41 (0.09) 


10.0 (4.2) 


35.7 (8.5) 


0.01 


88405.687 


H(52)/3 


1.01 (0.07) 


14.2 (0.9) 


27.2 (2.2) 


0.03 


91663.133 


H(64)(5 


0.52 (0.08) 


18.0(2.3) 


26.7 (4.5) 


0.02 


92034.434 


H(41)a 


4.01 (0.18) 


12.8 (0.7) 


33.9(1.8) 


0.11 


93607.316 


H(51)/3 


1.06(0.08) 


14.5 (1.3) 


34.2 (3.0) 


0.03 


93775.871 


H(58)r 


0.56 (0.09) 


15.6 (2.9) 


33.8 (5.8) 


0.016 


103914.838 


H(56)r 


0.46 (0.07) 


18.3 (2.7) 


40.1 (7.2) 


0.01 


105301.857 


H(49)/3 


0.66 (0.06) 


17.1 (1.3) 


27.3 (3.3) 


0.02 


106737.357 


H(39)a 


2.33 (0.07) 


16.4 (0.5) 


31.5(1.1) 


0.07 






MP2 (0' 


',40") 






Freq(MHz) 


Line 


I(Kkms-') 


v,,,(kms"') 


Av(kms-') 


Tmb(K) 


85688.390 


H(42)a 


0.89(0.16) 


11.0(2.1) 


23.2 (4.8) 


0.04 


88405.687 


H(52)/3 


0.47 (0.06) 


9.3 (2.2) 


35.0 (5.3) 


0.01 


92034.434 


H(41)a 


1.39(0.05) 


6.7 (0.6) 


33.4(1.4) 


0.04 


93607.316 


H(51)/3 


0.35 (0.08) 


6.6 (2.7) 


29.6 (10.7) 


0.01 


105301.857 


H(49)y8 


0.39 (0.05) 


9.4 (2.4) 


33.4 (4.4) 


0.01 


106737.357 


H(39)a 


0.57 (0.06) 


8.0(1.5) 


26.6 (3.2) 


0.02 


209998.777 


C(44)7 


0.51 (0.11) 


19.1 (0.9) 


8.2(1.8) 


0.06 



Note: " Overlapped with the HCO^ J=l— >0 line and self-absorbed; these two recombination lines are overlapped and the fit of the two each lines 
is uncertain; ' Doubtful detection; ^ This recombination line is contaminated by the emission of C"*0 J=2^1 from the image band. 
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Fig. A.l. Identified lines in offset (0",0")- In blue lines detected in signal band, in green lines identified in image band, and in red 
lines detected but not identified. 
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Fig. A.2. Identified lines in off'set (+15", -15")- In blue lines detected in signal band, in green lines identified in image band, and in 
red lines detected but not identified. 
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Fig. A.3. Identified lines in offset (0",40")- In blue lines detected in signal band, in green lines identified in image band, and in red 
lines detected but not identified. 
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Appendix B: Rotational Diagrams 
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Fig.B.l. Rotational diagrams of the IF's most representative molecules (Offset (0", 0"))- 
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Fig.B.2. Rotational diagrams of the MP's most representative molecules (Offset (+15", -15")). 
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Fig.B.3. Rotational diagrams of the most representative molecules towards MP2 (Offset (+0", +40"))- 



